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ABSTRACT 

It has been known for a long time that a large fraction of disc galaxies are lopsided. We 
simulate three different mechanisms that can induce lopsidedness: flyby interactions, 
gas accretion from cosmological filaments and ram pressure from the intergalactic 
medium. Comparing the morphologies, HI spectrum, kinematics and m = 1 Fourier 
components, we find that all of these mechanisms can induce lopsidedness in galax- 
ies, although in different degrees and with observable consequences. The timescale 
over which lopsidedness persists suggests that flybys can contribute to ~ 20 per cent 
of lopsided galaxies. We focus our detailed comparison on the case of NGC 891, a 
lopsided, edge-on galaxy with a nearby companion (UGC 1807). We find that the 
main properties of NGC 891 (morphology, HI spectrum, rotation curve, existence of 
a gaseous filament pointing towards UGC 1807) favour a fiyby event for the origin of 
lopsidedness in this galaxy. 

Key words: methods: A'^body simulations - galaxies: kinematics and dynamics - 
galaxies: interactions - galaxies: individual: NGC 891 



1 INTRODUCTION 

A high fraction of disc galaxies are lopsided, i.e. their gas 
and/or stellar component extend further out on one side 
of the galaxy than on the other (Baldwin, Lynden-Bell & 
Sancisi 1980; Block et al. 1994; Richter & Sancisi 1994; Rix 
& Zaritsky 1995; Schoenmakers, Franx & de Zeeuw 1997; 
Zaritsky & Rix 1997; Matthews, van Uriel & Gallagher 1998; 
Haynes et al. 1998; Swaters et al. 1999; Bournaud et al. 2005, 
hereafter B05; see Sancisi et al. 2008 for a review). 

The gaseous component of the disc is particularly af- 
fected by this phenomenon. Richter & Sancisi (1994) show 
that the lopsidedness of a galaxy can be inferred from asym- 
metries in its global HI profile, and estimate, from the anal- 
ysis of 1700 HI spectra, that ^ 50 per cent of disc galaxies 
are lopsided in the gaseous component. Haynes et al. (1998) 
confirm this result by the analysis of high signal-to-noise 
HI spectra of 104 galaxies, and suggest that some of the 
asymmetries may be induced by nearby companions (e.g. 
NGC 5324). Matthews et al. (1998) indicate that the inci- 
dence of gas lopsidedness is higher in the late-type galaxies 
(~ 77 per cent). The kinematics of the gas is often affected 
by lopsidedness: Swaters et al. (1999) find that the rotation 
curve of lopsided galaxies is rising more steeply on one side 
than on the other. 

Rix & Zaritsky (1995) and Zaritsky & Rix (1997), using 
near-infrared photometry of nearly face-on spiral galaxies, 
show that even the stellar component is lopsided in ~ 30 
per cent of their sample. Similarly, Rudnick & Rix (1998), 



using 7?— band photometry, find that ~ 20 per cent of their 
sample of nearly face-on early-type disc galaxies (SO to Sab) 
is lopsided in the stellar component. Thus, the incidence of 
stellar lopsidedness is similar for late-type and for early-type 
disc galaxies, although slightly lower in the latter case. The 
analysis of 25155 lopsided galaxies from the Sloan Uigital 
Sky Survey (Reichard et al. 2008) confirms that the lop- 
sided distribution of stellar light is due to a corresponding 
lopsidedness in the stellar mass. Finally, images and spa- 
tially integrated spectra of late-type galaxies (Rudnick, Rix 
& Kennicutt 2000) suggest a correlation between star for- 
mation and lopsidedness. 

The hypothesis that lopsidedness is due to galaxy in- 
teractions has been long discussed. Based on optical images, 
Odewahn (1994) finds that 71 of 75 lopsided Magellanic spi- 
rals have a nearby companion. However, Wilcots & Prescott 
(2004) obtain HI data of 13 galaxies from Odewahn (1994) 
and show that only four of them have Hl-detected neigh- 
bours. Thus, either lopsidedness is not related to galaxy 
interactions, or the asymmetries produced by these inter- 
actions are long-lived (surviving for ^ 2 orbital times after 
the encounter) and the lopsidedness persists even when the 
companion is quite far-off. 

From the theoretical point of view, the N-body simu- 
lations by Walker, Mihos & Hernquist (1996) suggest that 
minor mergers can induce lopsidedness over a long timescale 
( ^ 1 Gyr). However, B05 indicate that the lopsidedness 
produced by minor mergers disappears when the compan- 
ion is completely disrupted. Since most of observed lopsided 



2 Mapelli et al. 



galcixies are not undergoing mergers, the minor-merger sce- 
nario does not seem viable. BOS indicate that the most 
likely mechanism to produce lopsidedness is the accretion 
of gas from cosmological filaments. Alternative models sug- 
gest that baryonic lopsidedness can be induced by a lopsided 
dark matter halo (Jog 1997, 2002; Angiras et al. 2007) or by 
the fact that the disc is off-centre with respect to the dark 
matter halo (Levine & Sparke 1998; Noordermeer, Sparke & 
Levine 2001). 

In this paper, wc address the problem of the origin of 
lopsidedness by means of N-body/smooth particle hydrody- 
namics (SPH) simulations. In particular, we re-analyze in 
more detail the hypothesis of gas accretion, already pro- 
posed by B05, and we consider two new possible scenarios: 
the role of flyby interactions with smaller companions and 
that of ram pressure from the intergalactic medium (IGM) . 

For a comparison with observational data, we focus on 
the case of the edge-on galaxy NGC 891. We stress that 
quantifying lopsidedness in edge-on galaxies is more diffi- 
cult than in face-on galaxies, as bright regions on one side 
of the disc can be confused with lopsidedness. However, the 
lopsidedness of NGC 891 is well assessed (Sancisi & Allen 
1979; Baldwin et al. 1980; Rupen 1991; Swators, Sancisi & 
van der Hulst 1997). Furthermore, for NGC 891 recent HI 
observations are available, among the deepest ever obtained 
for an external galaxy (Oosterloo, Fraternali & Sancisi 2007, 
hereafter O07). This galaxy also shows many interesting pe- 
culiarities, e.g. the existence of a gaseous filament extending 
up to ~ 20 kpc vertically from the disc and located at ~ 10 
kpc from the centre of the galaxy. Finally, NGC 891 has 
also a smaller, gas-rich companion, UGC 1807, located at a 
projected distance of ~ 80 kpc, in the direction of the above 
mentioned gaseous filament. 



2 MODELS AND SIMULATIONS 

In this paper, we simulate three different processes: i) flyby 
interactions; ii) accretion from gaseous filaments; iii) ram 
pressure from the IGM. For all these scenarios we use a 
galaxy model similar to NGC 891 and whose main proper- 
ties are listed in Table 1. Such galaxy model has been gener- 
ated by using the method already described in Mapelli 2007 
(hereafter M07; see also Hornquist 1993; Mapelli, Ferrara & 
Rea 2006; MapelU et al. 2008a, 2008b). Here we briefly sum- 
marize the most important points, referring to M07 for the 
details. The galaxy model has four different components: 

- a Navarro, Frenk & White (1996, NFW) dark matter 
halo with virial mass M^ir = 1.4 X 10" Mq (O07), virial 
radius _R2()o = 104 kpc, and concentration c = 12; 

- a stellar exponential Hernquist disc (Hernquist 1993; 
M07) with mass Ma = 10^° Mq, scalelength Ra = 4.4 kpc 
(Shaw & Gilmore 1989), and scaleheight zo = 0.1 Rd.; 

- a stellar spherical Hernquist bulge (Hernquist 1993; 
M07) with mass M), = 2 x 10^ M©, and scalelength a = 
0.1 Ra; 

- a gaseous exponential Hernquist disc (Hernquist 1993; 
M07) with mass Mg = 4.1 x lO'' Mq (O07), scalelength Rg = 
Rd, and scaleheight Zg = 0.1 Rg. The gas is allowed to cool 
down to a temperature of 10* K, and to form stars according 
to the Schmidt law (Katz 1992). 



Table 1. Initial pajrameters for the model of NGC 891 and of 
UGC 1807. 





NGC 891 


UGC 1807 


M^ir/MQ 


1.4 X 10" 


7.65 X 10^ 


-R2()()/kpc 


104 


41 


Md/Mo 


10^° Mq 


9 X 10^ 


ild/kpc 


4.4 


0.5 


Zd 


0.1 Rd 


0.1 Rd 


MJMq 


2 X 10^ 





a 


0.1 Rd 





Mg/MQ 


4.1 X 10'^ 


4.5 X 10* 


Rg 


Rd 


Rd 




0.1 Rg 


0.1 Rg 



The model of NGC 891 has 619 500 halo dark matter 
particles, 500 000 stellar disc particles, 100 000 bulge parti- 
cles and 205 000 gaseous disc particles. Dark matter particles 
have mass equal to 2.23 x 10^ Mq, whereas disc, bulge and 
gas particles have mass equal to 2.23 x 10* M©. Softening 
lengths are 0.2 kpc for halo particles and 0.1 kpc for disc, 
bulge and gas particles (for the criterion used to estimate the 
softening, sec Dolmen 2001). The initial smoothing length of 
the gas is 0.1 kpc. 

The simulations have been carried out with the paral- 
lel A^-body/SFH code GASOLINE (August 2005 version; 
Wadsley, Stadel & Quinn 2004) on the cluster zbox2^ at the 
University of Ziirich. 

2.1 Flyby interaction 

In the case of flyby interactions, we simulate an intruder 
galaxy with the properties of UGC 1807 (see Table 1), i.e. 

- a NFW dark matter halo with virial mass M^ir = 
7.65 X lO'' A/0 (O07), virial radius R200 = 41 kpc, and con- 
centration c = 12; 

- a stellar exponential Hernquist disc with mass Md = 
9 X 10® Mq, scalelength Rd = 0.5 kpc, and scaleheight zo = 
0.1 Rd-, 

- a gaseous exponential Hornquist disc with mass Mg = 
4.5 X 10* M© (O07), scalelength Rg = Rd, and scaleheight 
Zg = 0.1 Rg. The recipes for gas cooling and star formation 
axe the same as for the target galaxy. 

As the mass of dark matter and baryonic particles in the 
intruder are 2.23 x 10^ M© and 2.23 x 10* M©, respectively 
(i.e. the same as in the target galaxy), the intruder is initially 
composed by 34 300 dark matter particles, 40 300 stellar disc 
particles and 20 200 gas particles. Also the softening lengths 
arc the same as in the target. 

We made various check runs, in order to find the initial 
centre-of-mass position and velocity of the intruder for which 
its final position and velocity best match the observations of 
UGC 1807. This best match is achieved for an initial centre- 
of mass position {x, y, 2;) = (— 38, 38, 80) kpc, and for an 
initial centre-of mass velocity {vx, Vy, Vz)={130, —130, 200) 
km s~^. 

^ http://www-theorie.physik.iiiiizh.ch/-dpotter/zbox/ 
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Table 2. Parameters of runs of gas accretion from filaments 



run 


coplanar 


2-shift" 


inclination'' 


ACCl 


Yes 








ACC2 


No 


10 kpc 





ACC3 


No 





-20° 



'^Off-set of the filament along the z-axis with respect to the 
plane of the galaxy. 

° Inclination of the filament with respect to the plane of the 
galaxy. 

2.2 Gas accretion from filaments 

The existence of filaments of cold gas that accrete onto 
galaxies is predicted by cosmological SPH simulations (Katz 
& White 1993; Katz et al. 1994; Keres et al. 2005; Ocvirk, 
Pichon & Teyssier 2008). It has also been proposed as a 
mechanism for the origin of polar ring galaxies (Maccio et 
al 2005). B05 showed that gas accretion might induce lopsid- 
edness, similar to the one observed in isolated spiral galaxies 
such as NGC 1637. 

We set up our simulations in a similar way to B05, i.e. 
we model a cylindrical, uniform gas filament which accretes 
onto the galax;y with a given accretion rate Mq ( ^ 6Mq 
yr~^). The filament is corotating with the galaxy. Different 
from B05, who consider only coplanar filaments, we also 
simulate the case of non-coplanar accretion. In the following, 
we will present the results for three different runs in which 
gas filaments are simulated (see Table 2). In all of them 
the filament has a radius of 12.5 kpc, a relative velocity of 
infalling onto the galaxy Vinfaii = 100 km s~^, an accretion 
rate Mq = 6Mq yr~^, and a 25 kpc off-set with respect to 
the centre of the galaxy (see fig. 16 of B05) . 

In run ACCl the filament is coplanar to the galaxy (see 
fig. 16 of B05). In run ACC2 the filament is still parallel 
to the plane of the galaxy, but is shifted of 10 kpc above it. 
Finally, in run ACC3 the filament has an inclination of —20° 
with respect to the plane of the galaxy. The existence of non- 
coplanar filaments, and in particular of the case ACC3, is 
supported by cosmological simulations (Keres et al. 2005). 

The total number of particles in the filament in the 
above runs is 543 910 and the mass of each particle is 2.23 x 
10'' Mq . We made also other check runs with smaller Mq 
(down to 2Mq yr~^). We integrate each run for ~ 1 Gyr. 

2.3 Ram pressure from the IGM 

The characteristic signatures of ram pressure have recently 
been observed also in dwarf irregular galaxies of the Local 
Group (McConnachie et al. 2007). This suggests the exis- 
tence of a non-negligible amount of diffuse IGM even in poor 
groups. In particular, gas is ram-pressure stripped from a 
galaxy if the density of the IGM is (Gunn & Gott 1972; 
McConnachie et al. 2007) 

>„„ ,„-6 -3 f 100 kms-^ V / Shi \2 
niGM^ 3.7x10 cm ^^^^ j ( 10"^ cm-^ J '^^^ 

where Vrei is the relative velocity between the galaxy and 
the IGM and Ehi is the column density of HI. 




Figure 1. Density map of gas, projected along the y-axis, for the 
target galaxy (edge-on) and the companion (face-on), at t = 640 
Myr after the beginning of the simulation and ~ 300 Myr after the 
fiyby. The frame measures 105 and 129 kpc, along the short and 
the long side, respectively. The density goes from 2.23 X 10"'^ 
(blue on the web) to 2.23 X 10^ Mq pc'^ (red on the web) in 
logarithmic scale. 

Thus, if the IGM density is higher than the threshold 
value in equation (1), even galaxies which are isolated or in 
small groups can suffer ram pressure and the distribution of 
their gas might appear lopsided. 

According to Bland-Hawthorn, Freeman & Quinn 
(1997), ram pressure heating can explain the fact that some 
galaxies (e.g. the lopsided galaxy NGC 253) have a truncated 
HI disc but present high temperature, ionized gas extending 
beyond the HI disc. Similarly, Ryder et al. (1997) suggest 
that ram pressure is responsible for the HI 'wake' observed 
in the spiral galaxy NGC 7421. 

To check the hypothesis that lopsidedness is induced 
by ram pressure, we made a simulation where the galaxy is 
moving into a cylinder of uniform gas with density mcM ~ 
5 X 10~^ cm~^ and Vrei ~ 200 km s~^. We made differ- 
ent runs, changing the inclination of the relative velocity 
Vrei with respect to the plane of the galaxy between and 
45°. The total number of IGM particles in the above runs 
is 2 088 630 and the mass of each particle is 2.23 x 10'' Mq. 
We integrate each run for ~ 1.5 Gyr. 



3 RESULTS 

The main result of this paper is that all the considered pro- 
cesses induce lopsidedness in the gaseous component. In the 
following, we discuss the details and the differences among 
the three scenarios. We firstly describe the morphological 
features of the simulated gaseous and (global and young) 
stellar component in the three scenarios (Sections 3.1-3.3). 
Then, we provide a more quantitative estimate of some ob- 
servational quantities, i.e. the Fourier m = 1 component 
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(Section 3.4), the HI spectrum and the rotation curve (Sec- 
tion 3.5). 

3.1 Flybys 

Fig. 1 shows the density of the gas in the target galaxy 
(seen edge-on) and in the companion galaxy (face-on) ~ 
300 Myr after the interaction, when the projected distance 
between the two galaxies is ^ 80 kpc, i.e. comparable with 
the observations of NGC 891 and UGC 1807. 

The simulation matches some of the properties of the 
NGC 891— UGC 1807 system. First, a gaseous filament con- 
nects the two galaxies in the simulation. The filament starts 
at ~ 10 kpc (projected distance along the y-axis) from the 
centre of the target galaxy, which is the position where the 
intruder reached the closest approach with the target. This 
filament is mostly due to the gas stripped from the compan- 
ion galaxy. The position of the filament is the same as in the 
observations (O07). However, in the simulation the density 
of gas is more or less constant along the filament, whereas 
in the data only the part of the filament which is closest to 
NGC 891 is visible. 

Second, the edge-on disc of the target galaxy is clearly 
lopsided (see the zoom in Fig. 2). In particular, the disc ex- 
tends further out on the side of the galaxy which is opposite 
with respect to the intruder. Indeed, the lopsidedness is in- 
duced by the flyby, as the perturbation originated exactly 
where the intruder penetrated the disc of the target. At the 
time shown by Figs. 1 and 2 (i.e. ~ 300 Myr after the flyby), 
the galactic disc has completed a half-rotation, and the per- 
turbation appears on the other side of the disc with respect 
to its initial position. 

Third, Fig. 3 shows that even the stellar component of 
the simulated galaxy has some degree of lopsidedness, al- 
though less evident than in the gas. The observations also 
show that the stellar component of NGC 891 is slightly lop- 
sided, especially the young population (van der Kruit & 
Searle 1981; Rand, Kulkarni & Hester 1990; Hoopes, Wal- 
terbos & Rand 1999; Kamphuis et al. 2007a; Kamphuis et 
al. 2007b). We also ran check simulations with different disc 
scalelength Rd (between 2.2 and 8.8 kpc), but we did not 
find significant differences from the point of view of lopsid- 
edness. Then, the ffyby scenario seems to reproduce quite 
weU the lopsidedness of NGC 891. 

Furthermore, the simulations suggest that the pertur- 
bation induced by flybys is long-lived ( 500 Myr) , explain- 
ing the origin of lopsidedness even in galaxies which have no 
longer interacting companions. This might be important, as 
many studies do not flnd any correlation between lopsided- 
ness and nearby companions (Rix & Zaritsky 1995; Zaritsky 
& Rix 1997; Wilcots & Prescott 2004; B05). 

3.2 Gas accretion 

Gas accretion from cosmological filaments is also able to 
produce lopsidedness, as suggested by B05. However, the 
kind of lopsidedness that might be induced by gas accre- 
tion is quite different from the one connected with flybys. 
This can be seen from the comparison between Fig. 2 and 
Fig. 4. Both figures show the density of gas in the edge- 
on (upper panels) and face-on (lower panels) view of the 




Figure 2. Density map of gas in the target galaxy, at t = 640 
Myr after the beginning of the simulation and ~ 300 Myr after the 
flyby. Top panel: the target galaxy is seen edge-on and the density 
is projected along the y-axis (this panel is a zoom of Fig. 1). 
Bottom panel: the target galaxy is seen face-on and the density is 
projected along the 2-axis. The frames measure both 80 kpc per 
edge. The density goes from 2.23 X 10"^ to 2.23 X 10^ Mq pc'^ 
in logarithmic scale. 



galaxy, but Fig. 2 refers to the flyby scenario and Fig. 4 to 
the gas accretion from a coplanar filament (ACCl). In the 
case of gas accretion, the lopsidedness appears much more 
pronounced, especially in the face-on view. Of course, this 
depends on the chosen accretion rate (M = Q Mq yr~^ in 
ACCl) and on the duration of the accretion phase (720 Myr 
in Fig. 4). However, even smaller accretion rates (down to 
M — 2Mq yr~^) and different timescales for the accretion 
(0.5—1.0 Gyr) produce approximately the same morphol- 
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Figure 3. Density map of stars, projected along the y-axis, in 
the target galaxy, a.t t = 640 Myr after the beginning of the 
simulation and ~ 300 Myr after the fiyby. Top panel: all stars 
of the target galaxy. The frame measures 80 kpc per edge and 
the density goes from 2.23 X lO'^ to 2.23 X 10^ Mq pc-^ in 
logarithmic scale. Bottom panel: young stars (< 100 Myr) of the 
target galaxy. The frame measures 80 and 22 kpc, along the long 
and the short side, respectively. The density goes from 7.05 X 10~® 
to 2.23 X 10^ Mq pc^^ in logarithmic scale. 



ogy. Therefore, the simulated galaxy looks more like, e.g., 
NGC 1637 than NGC 891. 

Of course, our models suffer from various limitations, 
such as the intrinsic problems of the SPH treatment of gas 
(see Agertz et al. 2007), the mass resolution and the fact 
that we do not account for various feedback processes (e.g. 
the cooling for temperatures below 10* K). Thus, we can- 
not take the morphology inferred from our simulations as 
a conclusive proof. However, the morphological difference 
seen in our simulations might suggest an intrinsic morpho- 
logical difference in observed galaxies between lopsidedness 
induced by gas accretion and lopsidedness induced by flybys. 
Thus, it would be interesting to compare the morphologies 
of observed lopsided galaxies, with or without companions, 
in order to check whether any difference can be found (as an 
example, we suggest the cases of NGC 1637 and NGC 891). 

Furthermore, the total stellar component (top panel of 
Fig. 5) does not seem to be lopsided in the gas accretion sce- 
nario, at odds with the observations of NGC 891. Instead, 
the young stellar population (i.e. the simulated star particles 
which have an age < 100 Myr) shows a clear lopsidedness 
(bottom panel of Fig. 5). This is likely due to the fact that 
gas accretion does not affect stars and induces lopsidedness 
only in the gas component. The distribution of young stars is 




Figure 4. Density map of gas in run ACCl, at t = 720 Myr after 
the beginning of the simulation. Top panel: the galaxy is seen 
edge-on and the density is projected along the j;-axis. Bottom 
panel: the galaxy is seen face-on and the density is projected 
along the z-axis. The frames measure both 100 kpc per edge. 
The density goes from 7.05 X 10"^ to 2.23 X 10^ Mq pc''^ in 
logarithmic scale. 

also lopsided, because these stars formed after the gas com- 
ponent has become lopsided. We can consider this process a 
sort of 'star formation' driven lopsidedness (see Section 3.4 
for details). 

The gas accretion from a coplanar filament does not 
produce anything similar to the gaseous filament observed 
in NGC 891. In runs ACC2 and ACC3 we change the orienta- 
tion of the filament (see Table 2), in order to see whether we 
can reproduce the filament observed in NGC 891. However, 
in runs ACC2 and ACC3 (Figs. 6 and 7, respectively) the 
galaxy disc becomes not only lopsided but also strongly dis- 
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Figure 5. Density map of stars in run ACCl, at i = 720 Myr 
after tlie beginning of tlie simulation. The galaxy is seen edge- 
on and the density is projected along the y-axis. Top panel: all 
stars. The frame measures 100 and 66 kpc, along the long and 
the short side, respectively. The density goes from 7.05 X lO'^ to 
2.23 X 10^ Mq pc"'^ in logarithmic scale. Bottom panel: young 
stars (< 100 Myr). The frame measures 100 and 27 kpc, along 
the long and the short side, respectively. The density goes from 
7.05 X 10~^ to 2.23 x 10^ Mq pc'^ in logarithmic scale. 



torted. This feature is not present in NGC 891 and in most 
of observed lopsided galaxies, whereas it might be somehow 
connected with warps. Such distortion tends to disappear 
only if the accretion rate is sufficiently small ( ^ 2 Mq yr~^). 
As cosmological simulations indicate that non-coplanar fil- 
aments exist (Keres et al. 2005), this may put limits on 
the density and/or on the lifetime of cosmological filaments. 
Further study of gas filaments in cosmological simulations 
is required, in order to quantify this effect. 

On the other hand, non-coplanar gas accretion does 
not produce any structure like the filament observed in 
NGC 891, unless we assume that the cosmological filament 
itself (or its remnant) is still observable today and represents 
the feature observed in NGC 891. Indeed, the co-existence 
of two different cosmological filaments, one coplanar induc- 
ing lopsidedness and the other non-coplanar producing the 
observed filament in NGC 891, might explain the morphol- 
ogy of this galaxy, but seems quite fine-tuning. 



3.3 Ram pressure 

Even ram pressure induces lopsidedness (see Fig. 8), after a 
sufficiently large amount of time ( ^ 600 — 700 Myr) . How- 
ever, lopsidedness produced by ram pressure appears less 
pronounced than that induced by either fiybys or gas ac- 
cretion. This is mainly due to the fact that ram pressure 
perturbs especially the low density gas in the peripheral re- 
gions, but does not affect the central parts of the disc. Thus, 
the effect of ram pressure is that of producing a faint tidal 
tail, more than a strong asymmetry. The stellar component. 




Figure 6. Density map of gas in run ACC2, at t = 720 Myr after 
the beginning of the simulation. Top and bottom panels are the 
same as in Fig. 4. The frames measure both 100 kpc per edge. 
The density scale is the same as in Fig. 4. 



both old and young, does not suffer significant perturbations 
due to ram pressure (see Fig. 9). 

Obviously, this result depends on our initial parameters, 
i.e. the density of the IGM and the relative velocity. On the 
other hand, a density higher than wigm = 5 x 10~^ cm"'' 
seems unrealistic for isolated galaxies (but this might be 
due to the lack of dedicated observations, as the paper by 
McConnachie et al. 2007 suggests) . As a lower density of the 
IGM would produce smaller perturbations, we infer that ram 
pressure cannot induce strong lopsidedness. 

Another point to stress is the role of the orientation of 
the relative velocity. We explored different inclinations of 
the relative velocity vector with respect to the plane of the 
galaxy, from to 45°. Large inclinations ( ^ 20°) produce 
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Figure 7. Density map of gas in run ACC3, at f = 720 Myr after 
the beginning of the simulation. Top and bottom panels are the 
same as in Fig. 4. The frames measure both 100 kpc per edge. 
The density scale is the same as in Fig. 4. 



tidal tails which are non-coplanar with the galaxy, and this 
does not look like lopsidedness. Only for smaller inclinations 
the tidal tail resembles a feature of lopsidedness. 

Finally, we did not find any way to reproduce the fila- 
ment observed in NGC 891 with ram pressure. Ram pressure 
with non-zero inclination with respect to the plane of the 
galaxy (e.g. 10°, as shown in Fig. 8) produces some vertical 
perturbations in the gas density, but these are transient and 
much smaller than the observed filament in NGC 891. 




Figure 8. Density map of gas in the ram pressure scenario, at 
t = 1 Gyr after the beginning of the simulation. Top panel: the 
galaxy is seen edge-on and the density is projected along the x- 
axis. Bottom panel: the galaxy is seen face-on and the density is 
projected along the z-axis. The frames measure both 100 kpc per 
edge. The density goes from 5.60 X 10~^ to 2.23 X 10^ Mq pc'^ 
in logarithmic scale. The two arrows indicate approximately the 
direction of the relative velocity between galaxy and IGM (which 
has an inclination of 10° with respect to the plane of the disc). 



3.4 m = 1 asymmetry 

In order to quantify our results, we derive the value of the 
Fourier m = 1 component for the surface density of our 
gaseous and stellar discs, m = 1 asymmetries are the best 
indicator of lopsidedness in stellar discs (Rix & Zaritsky 
1995; Zaritsky & Rix 1997; BOS). 

Fig. 10 shows the normalized strength of the Fourier 
component m = 1 (Al) as a function of radius, for all the 
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Figure 9. Density map of stars in the ram pressure scenario, at 
t = 1 Gyr after the beginning of the simulation. The galaxy is seen 
edge-on and the density is projected along the x-axis. Top panel: 
all stars. The frame measures 100 and 66 kpc, along the long and 
the short side, respectively. The density goes from 5.60 X 10"^ to 
2.23 X 10^ Mq pc~^ in logarithmic scale. Bottom panel: young 
stars (< 100 Myr). The frame measures 100 and 27 kpc, along 
the long and the short side, respectively. The density goes from 
7.05 X 10~^ to 2.23 x 10^ Mq pc'^ in logarithmic scale. 



considered models. Al has been defined as in equation (1) of 
BOS. We calculated Al for both the stellar and the gaseous 
component. In the case of the stellar component we distin- 
guish between the global stellar population (central panel of 
Fig. 10) and the young stars (bottom panel of Fig. 10). The 
behaviour of Al for the gaseous component is shown in the 
top panel of Fig. 10. In the case of gas accretion (ACCl, 
triangles in Fig. 10) Al rises very steeply, and reaches the 
value of~ 0.3 at 2 — 2.5 Rg. Al rises also in the flyby (circles 
in Fig. 10) and in the ram pressure (squares in Fig. 10) sce- 
narios, indicating lopsidedness, but it reaches smaller values 
(-0.1). 

We also calculated the value of (Al), i.e. the average 
value of Al between 1.5 and 2.5 disc scalelengths (see B05). 
According to B05, a galaxy is 'lopsided' when (Al) > 0.05. 
For the gaseous component, we found (Al) ~ 0.08, 0.28 and 
0.06 for the flyby, gas accretion and ram pressure scenario, 
respectively. Thus, all these three scenarios induce lopsided- 
ness in the gaseous component, but the asymmetry produced 
by gas accretion is much stronger than the one due to fly- 
by s and ram pressure. This confirms what we found from a 
qualitative analysis in the previous sections. 

The central panel of Fig. 10 shows the behaviour of Al 
for the global stellar component. In the flyby scenario the 
lopsidedness of the stellar population is similar to the one 
of the gaseous component: Al rises almost monotonously 
from — at the centre to 0.1 at 2.5 Rd- This is due to 
the fact that the flyby affects both the gaseous and the stel- 
lar component. Instead, in the gas accretion and in the ram 
pressure scenarios Al remains almost flat through the entire 
disc. Similarly, for the stellar component, (Al) ~ 0.08, 0.02 
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Figure 10. Normalized strength of the Fourier component m = 1 
as a function of radius, for gas (top panel), global stellar con- 
tent (central panel) and young stellar population (bottom panel). 
Open circles connected by solid line (red on the web): flyby sce- 
nario (at t = 300 Myr after the interaction); open triangles con- 
nected by dotted line (blue on the web): gas accretion (ACCl, at 
t = 720 Myr); open squares connected by dashed line (green on 
the web): ram pressure (at t = 1 Gyr). 



and 0.01 for the flyby, gas accretion and ram pressure sce- 
nario, respectively. This suggests that the stellar population 
is not lopsided in the case of gas accretion and ram pressure. 

This flnding might seem inconsistent with the results 
by B05, where gas accretion is shown to induce lopsidedness 
also in the stellar component. However, fig. 17 of B05 shows 
that (Al) < 0.05 for t < 0.7 - 0.8 Gyr, in agreement with 
our results. A value of (Al) ~ 0.2 is reached only at t ~ 2 
Gyr in the simulations presented by B05. This means that, 
in the gas accretion scenario, there is a time delay between 
the epoch in which lopsidedness appears in the gaseous dis- 
tribution and the epoch in which even the stellar compo- 
nent starts to be lopsided. The extent of this delay depends 
mainly on the efficiency of star formation. 

This interpretation is confirmed by the bottom panel of 
Fig. 10, which shows the behaviour of Al for the young stel- 
lar population (< 100 Myr). In the case of gas accretion the 
lopsidedness of the young stellar population is much higher 
than that of the total stellar population: (Al) is ~ 0.09 and 
~ 0.02 for the young stars and for the total stellar mass, re- 
spectively. The difference between the distribution of young 
stars and that of the total stellar component in the gas ac- 
cretion scenario is also evident, qualitatively, in Fig. 5 (see 
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Section 3.2). In the case of flybys, the young stars are lop- 
sided ((^1) ~ 0.07) approximately in the same amount as 
the global stellar content {{Al) ~ 0.08). In the case of ram 
pressure, the young stars are not lopsided ((^41) ~ 0.03), as 
already noted for the entire stellar population. 

This supports the idea that there are two possible forms 
of stellar lopsidedness: in the first case the lopsidedness of 
the stars originates at the same epoch and for the same 
mechanism as the lopsidedness of gas (e.g. for a flyby), in 
the second case the lopsidedness of the stellar component 
is a consequence of the (pre-existing) lopsidedness of the 
gas (e.g. for gas accretion). In the latter case, lopsidedness 
is present only (or mostly) in the young stellar population, 
unless the system had enough time to evolve after the pro- 
cess which induced lopsidedness in the gas. For example, 
after ~ 1 Gyr from the beginning of gas accretion, the dif- 
ference between spatial distribution of old and young stellar 
population is large, whereas no significant difference in the 
lopsidedness of young and old stars is expected, if the gas 
accretion from cold filaments started much more than ~ 2 
Gyr ago (and cosmological simulations suggest that cold gas 
accretion was more important at high redshift than today, 
Keres et al. 2005). 

Most of observational work on stellar lopsidedness (Rix 
& Zaritsky 1995; Zaritsky & Rix 1997; Rudnick & Rix 1998) 
focuses on quite 'red' bands (especially / and K), in order 
to track the mass distribution. As young stars account only 
for a small fraction of the light in these bands (especially in 
early- type disc galaxies), we expect that the observed stellar 
lopsidedness is different from what we find in the gas ac- 
cretion scenario, unless gas accretion started sufficiently far 
ago. Furthermore, observational studies indicate that lop- 
sidedness in the stellar light distribution is primarily tracing 
lopsidedness in the total stellar mass distribution (Rudnick 
& Rix 1998; Reichard et al. 2008). In particular, the inci- 
dence of stellar lopsidedness is similar in late-type (~ 30 
per cent, Rix & Zaritsky 1995) and in early-type disc galax- 
ies (~ 20 per cent, Rudnick & Rix 1998), indicating that 
the asymmetric light distribution is not due to recent asym- 
metric star formation, but refiects global asymmetric mass 
distribution. This implies that the observed stellar lopsided- 
ness may be connected with fiyby events or with long-lived 
gas accretion, but not with recent ( 1 Gyr) gas accretion. 

On the other hand, Zaritsky & Rix (1997) and Rudnick 
et al. (2000) show that there is a correlation between lop- 
sidedness and star formation rate. This suggests that, even 
if lopsidedness is not due to recent star formation, the event 
which produced lopsidedness has also triggered star forma- 
tion. This matches quite well both the fiyby and the gas 
accretion scenario, where a star formation enhancement by 
a factor of ;i 2 is induced by the galaxy interaction and by 
the feeding of fresh cold gas, respectively. 

3.5 Spectrum and kinematics of gas 

As we discussed in the previous section, all the considered 
mechanisms induce a certain degree of lopsidedness in the 
density distribution of disc galaxies. We now compare the 
effects that these three mechanisms have on the HI spectrum 
and on the kinematics, which in the observational data also 
appear to be affected by lopsidedness. 

In order to derive the HI spectrum and the rotation 
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Figure 11. HI spectra of the simulated galaxies compared with 
the observations of NGC 891. Open circles (green on the web): 
data from Richter & Sancisi (1994). Solid line (blue on the web): 
simulations. Top panel: fiyby interaction; central panel: gas ac- 
cretion (ACCl); bottom panel: ram pressure (case reported in 
Fig. 8). 



curve of the galaxy, we adopt a procedure as similar as 
possible to the one used by observers: we rotate the sim- 
ulated galaxy by the same inclination angle as the observed 
NGC 891 (i.e. ~ 89° ) and we calculate the velocity of gas 
particles along the line-of-sight (the details of this procedure 
are given in Mapelli et al. 2008b). 

The derived HI spectra are shown in Fig. 11 and com- 
pared with the observational data of NGC 891 (Richter & 
Sancisi 1994). These results are quite interesting, as the 
gas accretion scenario (central panel) has some problem in 
reproducing the data, whereas the ram pressure (bottom 
panel) and especially the fiyby scenario (top panel) match 
quite well the observations. This might be crucial, if not for 
all the lopsided galaxies, at least for the case of NGC 891. 
In fact, other lopsided galaxies have significantly different 
spectra (see fig. 3 of Richter & Sancisi 1994), which might 
be reproduced even by gas accretion (e.g. NGC 5006 or 
NGC 4303). However, the spectrum of NGC 891 seems to 
agree only with the scenarios of ram pressure and fiyby. 

Analogous conclusions can be inferred from the analysis 
of the rotation curves (Fig. 12). We calculated two 'partial' 
rotation curves for each simulated galaxy, one for the reced- 
ing side of the galaxy (circles in Fig. 12) and the other for 
the approaching side (squares in Fig. 12). This procedure 
is analogous to the position-velocity diagram and has been 
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4 DISCUSSION AND CONCLUSIONS 
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Figure 12. Rotation curves of the simulated galaxies for the 
approaching (open squares connected by dashed line, red on the 
web) and receding (open circles connected by solid line, blue on 
the web) sides separately. Top panel: fiyby interaction; central 
panel: gas accretion (ACCl); bottom panel: ram pressure (case 
reported in Fig. 8). 



adopted to quantify kinematic lopsidedness in Swaters et 
al. (1999). 

In the flyby scenario (top panel) the rotation curve de- 
clines gently at larger radii than ~ 5 kpc, especially in the 
receding side. The same trend is apparent from the new data 
of NGC 891 (see fig. 5 of O07). A similar behaviour can be 
found in the ram pressure case (bottom panel), although 
the decrease starts at larger radii and is less pronounced. 
Instead, in the gas accretion case (central panel) the rota- 
tion curve is almost flat everywhere. 

Furthermore, in the flyby scenario the rotation curve of 
the receding side is more steeply rising than the one of the 
approaching side in the range between ~ 5 and ~ 10 kpc. 
An analogous trend also exists in the partial rotation curves 
of NGC 4395 and DDO 9 (Swaters et al. 1999), while the 
case of NGC 891 is uncertain (O07). Also in the gas ac- 
cretion scenario the rotation curve of the receding side is 
steeper, but at larger radii ( ^ 6 — 7 kpc) , whereas in the 
ram pressure case the two sides behave in the same way. We 
stress that in some observed (Rubin, Hunter & Ford 1991; 
Hernandez-Toledo et al. 2003) and simulated (Pedrosa et al. 
2008) interacting galaxies the receding side of the rotation 
curve is slightly different from the approaching side. This 
circumstance may also suggest a connection between lopsid- 
edness and fly by s. 



In this paper we considered three different physical mech- 
anisms acting on disc galaxies: flyby interactions, accretion 
from cosmological gas fllaments and ram pressure from the 
IGM. We showed that all of these processes might induce 
lopsidedness in the gaseous component of discs. Likely, these 
three scenarios are not exclusive, but might account for lop- 
sidedness in different galaxies. 

Our simulations suggest that the features of lopsided- 
ness induced by these three processes are slightly differ- 
ent from each other. In particular, ram pressure creates 
only moderate tidal tails in galaxies, whereas flybys account 
for much stronger asymmetries and gas accretion from fll- 
aments might produce even more pronounced lopsidedness. 
The spectra and the kinematic features connected with these 
three scenarios are also slightly different. Furthermore, ram 
pressure is not able to induce lopsidedness in the stellar pop- 
ulation, whereas flybys simultaneously produce asymmetries 
in the gaseous and in the stellar disc. Gas accretion from fll- 
aments generates asymmetries only in the gaseous disc, but 
the star formation in the asymmetric gaseous disc induces 
lopsidedness in the stellar component at later times. Thus, 
our simulations suggest that there is not a unique manifesta- 
tion of lopsidedness but many different forms: lopsidedness 
may be driven by various mechanisms, and, depending on 
this, it may affect only the gas and not (or not yet) the stel- 
lar population, it may be more or less pronounced, it may 
have different timescales. This might also explain why the 
observed incidence of lopsidedness is different in the stel- 
lar and in the gaseous component. Thus, it is interesting 
to look at the observational properties of lopsided galaxies, 
and check whether their morphological features allow to in- 
fer which process produced their lopsidedness. In particular, 
the characteristics of mildly-lopsided (e.g. NGC 891) and 
strongly-lopsided (e.g. NGC 1637) galaxies can be studied 
separately, in order to flnd any intrinsic difference between 
them (e.g. differences in the incidence of companion galax- 
ies, in the star formation rate, in the distribution of young 
stars). Furthermore, lopsided galaxies are also observed in 
groups and clusters (Haynes, Giovanelli & Kent 2007, and 
references therein) where flyby interactions and ram pres- 
sure effects can be even more important. It would be in- 
teresting to compare the fraction of lopsided galaxies as a 
function of environment. 

Another possible feature of lopsidedness which needs to 
be investigated, both by theoretical models and by observa- 
tions, are metallicity gradients. There are no specific obser- 
vations of metallicity gradients in lopsided galaxies. In the 
case of both flybys and ram pressure no signiflcant metal- 
licity gradients are expected. Instead, in the case of cold ac- 
cretion the metallicity of the gaseous fllament is likely lower 
than the one of the target galaxy. For example, Ocvirk et 
al. (2008) derive a metallicity Z ~ 10"*^ - 10"^ Zq (with a 
very large spread) for the cold gas accreted from low density 
filaments at redshift z = A, whereas the median metallic- 
ity for a Milky Way-like galaxy is ~ 0.5 Z©. However, it is 
unclear whether this feature is preserved (and observable) 
in the nearby Universe. In fact, it is likely that stars form 
from a mixture between the gas in the cold filament and the 
gas which was already in the disc. Furthermore, if the gas 
accretion from the filament is not recent, but occurred in 
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the past, pollution from supornovae, stellar winds and other 
feedback rnochanisrns have already erased any metallicity 
gradient both in the gaseous and in the stellar component. 

Here, we focused our analysis on the CEise of NGC 891. 
Its morphology, the existence of a companion and of a 
gaseous filament pointing towards it, the observed HI spec- 
trum and also the shape of the rotation curve favour the 
hypothesis of a fiyby as source of lopsidedness. However, 
we cannot exclude other scenarios. In particular, the ori- 
gin of the gaseous filament observed in NGC 891 might bo 
completely different from that of lopsidedness (e.g. a galac- 
tic fountain; but see Praternali & Binney 2006). Other ob- 
served galaxies (e.g. NGC 2985, Sancisi ct al. 2008) show 
features of lopsidedness similar to those of NGC 891 and 
have nearby companions. For example, NGC 3027, the com- 
panion of NGC 2985, has an asymmetric HI distribution. 

Previous studies (Wilcots & Prescott 2004; B05 and ref- 
erences therein) tend to exclude a link between galaxy in- 
teractions and lopsidedness. However, B05 do not consider 
flybys but only minor mergers. The main difference between 
a merger and a flyby is that the merger remnant remains 
visible inside the target galaxy for a long time after the 
beginning of the merger. Thus, if lopsidedness was due to 
mergers, we should also detect the remnants, which is not 
the case^. Instead, in a flyby, the companion goes away after 
the interaction. As the effects of the flyby are visible for at 
least ~ 500 Myr after the interaction, the companion has 
enough time to travel far from the target galaxy. 

Previous observations do not find a significant correla- 
tion between lopsidedness and nearby companions. However, 
observations such as those by Wilcots & Prescott (2004) con- 
sider only companions which are still interacting with the 
target galaxies. In such analysis, companions which had a 
non-recent flyby with the target are ignored. Even NGC 891 
would be considered an isolated galaxy according to this 
method. Similarly, B05 measure the incidence of compan- 
ions by the tidal parameter (see equation 4 of B05). The 
tidal parameter of NGC 891 is ~ — 5.0, which classifies this 
galaxy as rather isolated (see fig. 9 of B05). 

Furthermore, most of the previous studies assume that 
lopsidedness is produced by the same process in all the 
observed galaxies. In this paper we showed that at least 
three different processes may be responsible for lopsided- 
ness. Thus, the problem of lifetime of lopsidedness induced 
by flybys disappears: it may be that long-lived lopsidedness 
is connected with gas accretion, and in this case the lopsided 
galaxy appears completely isolated; whereas short-lived lop- 
sidedness is connected with fiyby, and in this case the lop- 
sided galaxy has a companion, although non-necessarily an 
interacting companion. Thus, wo argue that the role of fly- 
bys in producing lopsidedness cannot bo excluded on the 
basis of current observations and simulations. 

Therefore, an important issue to address is the rate 
of flybys which may induce lopsidedness. Here we try to 
do an approximate calculation. Another class of galaxies is 
known to be produced by flybys: the coUisional ring galax- 
ies (CRGs, see Mapelli et al. 2008a, 2008b and references 



therein). The density of CRGs in the local Universe is 
ncRG ~ 5.4 X 10~® Mpc~^ (where h is the Hubble pa- 
rameter. Few & Madore 1986). Ring galax:ies are produced 
during interactions where the mass of the intruder is > 1/20 
of the mass of the target (Appleton & Struck-MarccU 1996, 
and references therein). From check simulations we found 
that lopsidedness is induced during galaxy interactions, if 
the intruder mass is at least 1/20 of the target mass. Thus, 
the mass-ratio between intruder and target is similar (within 
a factor of 2) for interactions producing ring galaxies and for 
encounters inducing lopsidedness. However, other character- 
istics of flybys which induce lopsidedness are different from 
those of interactions which produce CRGs. In particular, to 
produce a ring, the flyby must have a very small impact 
parameter (&crg ~ 0.36 -Rj, corresponding to ~ 8 per cent 
of the galaxy radius. Few & Madore 1986). Instead, lopsid- 
edness is produced by off-centre interactions, with impact 
parameter as large as 6max ~ 4i?d. Furthermore, regular 
CRGs cannot form when the inclination angle between the 
disc axis of the target and the intruder is too large (~ 45°, 
Few & Madore 1986), whereas flybys producing lopsidedness 
do not have this limitation. Finally, ring galaxies are short- 
lived (tcRG ~ 500 Myr, Mapelli et al. 2008b), whereas in our 
simulations the lopsidedness is still present at tiop ^ 1 Gyr 
after the flyby. Thus, we can derive the density of galaxies 
which are lopsided as a consequence of flybys {nay) by con- 
sidering the density of CRGs and by correcting it for these 
differences^ : 



rifly ~ ncRG ■ 



"CRG 



i'CRG 1 ftop 

1 — COS 45° tCRG 



(2) 



Substituting in equation (2) the values reported above, we 
derive nay ^ 4.5 x 10~^ Mpc~^. Since the density of rel- 
atively bright disc galaxies (i.e. with absolute magnitude 
;S -19 mag) is naisc ~ 4.9 x 10"^ Mpc"^ (Few & Madore 
1986 and references therein), this suggests that only one 
galaxy every ~ 11 spiral galaxies can bo lopsided as a con- 
sequence of flybys. As lopsidedness of the gas component is 
observed in ~ 50 per cent of disc galaxies (Richter & Sancisi 
1994), flybys can reasonably contribute to ~ 20 per cent of 
lopsided galaxies. 

Unfortunately, the density of lopsided galaxies which 
have undergone gas accretion or ram pressure is much more 
difflcult to estimate. For ram pressure, the peculiar veloci- 
ties of ~ 149 galaxies in the Local Volume (i.e. ^ 10 Mpc 
from the Milky Way) are known (Whiting 2005) and ^ 21 
per cent of them have velocity higher than 100 km s~^. How- 
ever, no data are available for the local density of the IGM, 
and we cannot even derive an upper limit for the density of 
ram-pressure perturbed galaxies. For gas accretion, cosmo- 
logical simulations show that most of galaxies accrete from 
cosmological filaments. Thus, all galaxies may have under- 
gone cold accretion. Those galaxies in which the accretion 
was predominantly coplanar may have become lopsided. 

We stress that accretion from non-coplanar filaments 
might produce distortions in the disc, different from lopsid- 
edness. This fact deserves further investigations, as it might 
put constraints on the lifetime, the accretion rate and the 



^ However, some lopsided galaxies (e.g. MlOl) show the signa- 
tures of mergers occurred in the past (see Sancisi et al. 2008 and 
references therein). 



^ We remind that the geometric cross-section of intera<;tions with 
impact parameter smaller than b is proportional to 
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existence itself of non-coplanar filaments. Furthermore, the 
existence of non-coplanar filaments might be associated with 
the formation of warps (Binney 1992). 

Finally, in this paper, we did not consider two further 
possible mechanisms for the origin of lopsidcdncss: the de- 
velopment of asymmetries due to lopsided dark matter halos 
(Jog 1997, 2002) or to displacements between disc and halo 
(Levine & Sparke 1998). These mechanisms might coexist 
with the three scenarios analyzed in this paper. However, 
these mechanisms cannot produce the filament observed in 
NGC 891. 
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